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Abstract. We prove that the level sets of a real C s function of two 
variables near a non-degenerate critical point are of class c' S//2 ' and apply 
this to the study of planar sections of surfaces close to the singular section 
by the tangent plane at an elliptic or hyperbolic point, and in particular 
at an umbilic point. 

We also analyse the cases coming from degenerate critical points, cor- 
responding to cusps of Gauss on surfaces, where the differentiability is 
now reduced to C' s/ ^ 4 '. 

However in all our applications to symmetry sets of planar sections of 
surfaces, we assume the C°° smoothness. 



1 Introduction 

The medial axis or skeleton of a closed plane curve 7 encode a great deal of in- 
formation about the shape of the region enclosed by the curve. This information 
has been exploited in many ways 1 . The medial axis can be defined, for a smooth 
curve 7, as the closure of the locus of centres of 'bitangent' circles — tangent to 7 
in two places — and whose radius equals the minimum distance from the centre to 
7. The symmetry set is the closure of the locus of centres of all bitangent circles. 
It thus has a richer structure which underlies the transitions of the medial axis. 
See for example 6 . 

The local structure of the symmetry set of a generic plane curve is of four 
kinds: smooth branches, endpoints, cusps and triple crossings. All of these can 
be seen in the example of Figure |21 (For medial axes there are only smooth 
branches, endpoints and Y-junctions.) Symmetry sets and medial axes of curves 
which vary in a generic 1-parameter family 7^ are well understood provided 
the curves of the family remain nonsingular; a complete local classification was 
given in [2}. In this paper, we consider instead the family of plane sections jk 
of a smooth surface M in 3-space, where the plane moves parallel to itself and 
becomes the tangent plane to M at a point p for say k = 0. Then the section 
70 will be singular at the tangency point p. For p an elliptic point 70 is, locally, 



1 Typing 'medial axis' into an internet search engine produces hundreds of references, 
to character recognition, Voronoi diagrams, interrogation, reconstruction, modifica- 
tion and design of shape, etc. 
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a single point. For p a hyperbolic point 70 is locally two smooth transversally 
intersecting branches, while at an ordinary parabolic point 70 has an ordinary 
cusp at p. See Figure ^ which also shows a 'hyperbolic cusp of Gauss' where 
the intersection is two smooth tangential branches. See H4.2I In this situation, 
the results of [2] are invalid and we need to use other techniques to find out how 
the symmetry sets (and medial axes) behave in the family of curves. 




Fig. 1. Surfaces with the tangent planes at elliptic, hyperbolic, parabolic and 
hyperbolic cusp of Gauss points. The intersections are a point, two smooth trans- 
verse curves, a cusped curve and two tangential curves respectively. Figure pro- 
duced with MAPLE. 

The motivation for this work comes from the study of isophote curves of a 
2D image, which are sections of an intensity surface in 3-space. Results on the 
pattern of vertices (extrema of curvature) and inflexions (zeros of curvature) of 
the sections 7& have been given in 0] while in [S] there are results on the patterns 
of cusps and triple crossings of the symmetry set. Let the surface M be given 
locally by an equation z = f(x,y). Then the sections jk are by their nature given 
by equations f(x,y) = k, and not by parametrizations j(t) = (X(t),Y(t)). Now 
plotting symmetry sets of parametrized curves is reasonably straightforward (see 
2J. But curves given by equations are a different matter altogether: there is not 
even any foolproof way to determine how many real components a curve f{x,y) = 
k has. Furthermore it may not be feasible to find an exact parametrization of 
one component of a curve f(x, y) = k. 

In this paper we give a result on parametrizing plane sections of surfaces and 
use it to approximate such a section to any desired degree of accuracy. We apply 
this to the plotting of symmetry sets and medial axes of the surface sections 7&. 
In [J21we give the main theoretical result and its proof (Propositions , in 
we give the way in which this is implemented in practice, and in jJHwe apply the 
method to symmetry sets, concentrating in ij4.ll on the case of umbilic points 
(elliptic points where the principal curvatures are equal). In i)4.2l we consider 
the case of an 'elliptic cusp of Gauss', for which the curves 7^. are simple closed 
curves when k > 0. 
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2 Level sets of functions 

The intuitive idea here that we 'blow up' the origin to turn the surface z = f(x, y) 
into one with nonsingular sections. We shall take the surface M in Monge form, 
that is with 

f(x, y) = \ (kix 2 + n 2 y 2 ) + b x 3 + b 1 x 2 y + b 2 xy 2 + b 3 y 3 
+ cqx 4 + cix 3 y + c 2 x 2 y 2 + c 3 xy 3 + c^y 4 

+ d x 5 + d 1 x i y + d 2 x 3 y 2 + d 3 x 2 y 3 + d 4 xy 4 + d 5 y 5 + h.o.t. (1) 

For simplicity, consider the case where k,\ and k 2 are > so that, locally to 
the origin, z > at points of M. Let us write x = tX, y = tY, z = t 2 ; after 
cancellation of t 2 the equation z — f(x,y) then becomes 

1 = i(«iX 2 + k 2 Y 2 ) + t(b X 3 + ...)+ < 2 (c X 4 + ...) + ..., 

which is a smooth surface in (X, Y, i)-space whose sections t — constant / are 
scaled versions of the sections of M. We proceed to give the formal details of 
this procedure in a more general setting and in the next section show how in 
practice we have used it. 

We prove that the level sets of a real function of two variables near a non- 
degenerate critical point are of class C^/ 2 ^ 1 if the function is of class C s . To 
this end we consider / 6 C S (D), Del 2 the open unit disc, /(0) = df(0) = 0, 
and set q — Vdf (0) to be the hessian of / at the origin. We assume that q is 
cither positive definite or indefinite. See Proposition 01 below for the case of a 
degenerate critical point. 

Proposition 1. Assume that q is positive definite and let r$ : S 1 — > R + , where 
S 1 C K 2 is the unit circle, be defined by r (9) = [<7(cos#,sin0)] -1 / 2 , 6 E S 1 
corresponding to e l ° . 

Then there exists a function r : [0, e) x S 1 — > M + for some e > with 

(a) r e C[ s / 2 l([0,e) x S 1 ), 

(b) r(O,0)=r„(0) onS\ 

(c) f(t r{t, 9) cos 6, t r{t, 9) sin 9) = t 2 on [0, ejxS 1 . 

Remark: The equation (c) implies that for fixed t £ [0, e), the parametrized 
curve 9 ^ (r(t, 9) cos 9, r(t, 9) sin 9) is a rescaled level curve of the function /. 
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Proof. For t > consider the curve 

C t = {f(tx,ty)=t 2 , z = t 2 }cR 3 xyz 
Consider the surface defined by 

Q = {f(VW\ x > VW\y) = \ z \ for and ifav) = 1 for z = °}- 

We claim that Q is a regular surface of class C^ 2 ' near {z — 0} flQ- Given this, 
we have Ct — Q H {z — t 2 }, and since the intersection is transversal (as is seen 
by computing grad f(y\z\ x, y/\z~\ y) — \z\ ), we conclude that C t £ C*[ s / 2 ] with 
uniform bounds. Using the implicit function theorem we derive the existence of 
r and hence (a) to (c). Note that now Ct is parametrized by 

{(< r{t,9)cos9,t r(t, 9) sin 9, t 2 ) for 6> £ S 1 }. 

□ 

The claim above follows from the following proposition: 

Proposition 2. Let h(x) be a real function of one real variable, h = 0(\x\ 2 ) 
and h G C l . Then g(x,y) = h(^/\y\ x) is of class C^ 2 1 and g = 0{\y\). 

Proof. The proof of this proposition follows from Taylor's formula. □ 
Finally, we have the analogous result for q indefinite, which is proved in a 
similar way: 

Proposition 3. Assume that q is indefinite and of the form q{x,y) = n\x 2 — 
K2y 2 with Ki > K2 > 0. Let S + C S 1 be defined by q\s+ > 0. For any So C 
S + , Sq S + there exists e > and r : [0, e) x So —> K + with 

(a) re^/ 2 l([0,e)x5 ), 

(b) / (t-±^r{t,9) cosh 0,t^=r {9, t) sank 6) = t 2 on S x [0,e) 

(c) r(6»,0) = 1 on S . □ 

We have a result similar to Proposition ^ in the case of an elliptic cusp 
of Gauss. In this case the quadratic terms of / are degenerate, equal to x 2 
say, so that the surface z = f(x,y) has a parabolic point at the origin. By 
making x divide the cubic terms we can still have a closed curve of intersection 
f(x, y) — k > provided the terms in x 2 , xy 2 , y give a positive definite quadratic 
form in x and y 2 . This is called an elliptic cusp of Gauss. The contact of the 
surface z = f(x,y) with its tangent plane at the origin is of type A3 in the 
notation of Arnold. See for example |1I8| for many geometrical properties of 
these points. 

Proposition 4. Assume that f £ C S (D) and let the point £ M be an elliptic 
cusp of Gauss and Ki > 0. Set q(x, y) = \k\x 2 + b2xy 2 + C4y 4 , with b\ < 2k,\C\. 
Let ro : S 1 — > ]R + be defined by ro(0) — (/(cos 9, sin 9)~ 1 / 4 . Then there exists a 
function r : [0, e) x S 1 — > R + for some e > with 



Level sets 5 



(a) r e Cl s / 4 ]([0,e) x S 1 ), 

(b) r(O,0)=r o (0) onS 1 , 

(c) f{{tr{t,9)) 2 cos9 1 tr{t,9)sm9) = t 4 on [0, e) x 5 1 . 

Proof. The rescaling here is nonhomogeneous and given as follows. Let 

C t = {/(i 2 My)=i 4 , ^i 4 KIR^ 2 . 

Consider the surface defined by 

Q = a;, |z| 3 y) = |;z| for z ^ and q(x,y) = 1 for # = 0}. 

Now Q is a regular surface of class C^'^.This follows from a result analogous 
to Proposition 2. □ 
Remark: For a fixed value of t, we are parametrizing the level set f(x, y) = t 4 
not in 'polar coordinates' where each ray from the origin intersects the curve in 
one point, but by means of parabolas of the form y 2 — kx for constants k. 



3 Finding the parametrization in practice 

All functions and surfaces from now on will be assumed smooth of class C°° , 
that is s — oo in the Propositions of SJU 

We seek to approximate the sections f(x,y) — k of the surface M up to a 
suitable order. Let M be given in Monge form JIJ. As in S]21the two cases (i) 
elliptic: Ki > 0, n 2 > and (ii) hyperbolic: K\ > 0, k 2 < are different: in the 
former case the sections f(x,y) = k are, locally to the origin, closed curves for 
small k > whereas in the latter case they are open curves extending to infinity 
for both signs of k. 

In this article our applications will concentrate on the cases where the inter- 
section is a closed curve, but we give some details of other cases in this section. 

• For a hyperbolic point we seek to parametrize the section by 

x = t r(t, 8) cosh^, y = t r(t, 9) sinh 9, z = t 2 , or 

x = tr(t,9)smh0, y = t r(t, 6) coshfl, z = -t 2 , (2) 

for a suitable function r. 

• For an elliptic point we seek to parametrize the section by 

x = t r(t, 6) cos 9, y = tr(t, 9) sin 9, z = t 2 , t > 0, < 9 < 2tt. (3) 

• For an 'elliptic cusp of Gauss' we seek to parametrize the section by 

x = (t r(t,9)) 2 cos9, y = tr(t,6)8ia6, z = t 4 , t > 0, < 9 < 2tt. (4) 



We will write 

r(t,9) ^r Q +r 1 t + r 2 t 2 + 
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where the r, are functions of 9 only. In the elliptic case for small values of k > 
the section z = k is, locally to the origin, a closed curve and the r, will be 
periodic and hence functions of cos# and sin#. 

We therefore need to express the in terms of the Monge form coefficients 
Ki, bi, Ci and so on. This is done by comparison of series, and enables us to 
approximate the surface z = f(x,y), and the sections z = constant up to any 
desired accuracy. 

Consider the terms of degree i in the Taylor expansion of /: this is a homoge- 
neous polynomial of degree i in x and y. Denote by pi the function of 9 obtained 
from this homogeneous polynomial as follows: 

hyperbolic case: replace x by cosh# and y by sinh# when z > and 

x by sinh# and y by cosh# when z < 0; 

elliptic case: replace x by cos# and y by sin#. Thus 

P2 = \(k x c 2 + k 2 s 2 ); 

p 3 = b a c 3 + b X (?s + b 2 cs 2 + b 3 s 3 ; 

P4 = c c 4 + ac 3 s + c 2 c 2 s 2 + c 3 cs 3 + C4C 4 ; 

p 5 = d c 5 + dic^s + d 2 c 3 s 2 + d 3 c 2 s 3 + d A cs 4 + d 5 s 5 ; 

and so on, 

where 

hyperbolic case: c — cosh#, s = smh9 if z > and c = sinh(9, s — cosh# if z < 0, 
elliptic case: c = cos 9 and s = sin#. 

Hyperbolic case By substitution in the Monge form 0J, we obtain in succes- 
sion the following formulas. 

rn = — \= when z > and rn = , 1 when z < ; 
r2 = ~2^( r oPi + 3rgrip 3 + r\p 2 )- 

r 3 = ~2F^( 2r i r 2P2 + r%p 5 + 3r^r 2 p 3 + 3r r 2 p 3 + Ar^r^); 

r4 = ^2^( 5r o r iP5 + ^r 2 p A + 2r 1 r 3 p 2 + §r\r\p± + r\p 2 + 6r a rxr 2 p 3 

+ 3r 2 r 3 p 3 + r\p 3 + r%p 6 ); 
r 5 = ^27^( 3r o r 2P3 + 3r 2 r 2 p 3 + 6r rir 3 p 3 + 3r 2 r 4 p 3 + 12r 2 rir 2 p 4 

+ 4r$r 3 p 4 + 4r rfp 4 + 6r^rip 6 + 5r 4 r 2 p 5 + 10r$r 2 p 5 + 2r 1 r A p 2 + 2r 2 r 3 p 2 ); 

etc. 

Taking the quadratic terms of the surface to be x 2 — a 2 y 2 where a > 0, we 
have 

rl = 3_ ( z > 0): rl = ^— (z < 0). 

(1 -a 2 ) cosh 2 + a 2 ° 1 - (1 -a 2 )cosh 2 <9 V ' 

When z > 0, the expression for 7q is > for all 9 when < a < 1, but for 

2 2 

a > 1 we need cosh 9 < " , , that is 1 < cosh# < . g . 

a — 1 y a- — 1 

When z < 0, the expression for r 2 . is > for all 9 when a > 1 but if < a < 1 
we need 1 < cosh 9 < . 1 , . 

VI — a 2 

Elliptic case This is much simpler because, as above, we can expect a global 
parametrization of the closed curves f(x, y) = k, locally to the origin, for k > 0. 
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The formulae for pi are given above and those for n are exactly the same as for 
the hyperbolic case, except that, for z > only, 



Ati cos 2 6 + n 2 sin 9 

Note that we are assuming Ki > 0, k 2 > since the surface is assumed to be 
above the plane z — close to the origin. Thus tq is always real. 

When the origin is an umbilic point, then k% = «2 = K , say, and r 2 , = -, a 
constant. Indeed, we always scale the variables so that the quadratic terms of / 
are x 2 + y 2 and then r 2 has the constant value 1 and we can take r = 1 without 
loss of generality. 

Elliptic cusp of Gauss case The relevant parametrization here is J2J, and we 
apply Proposition 0] For a cusp of Gauss, after renaming axes if necessary, we 
have 

, . j/ & 2< 2k i c 4 elliptic cusp, 

Kl * 0, n 2 = 0, b 3 =0and| 6i> ^ hyperbolic cugp (5) 

Consider the elliptic case. The expansion of the function r in Proposition 01 is 
obtained by a 'weighted' version of the method used in the previous cases. This 
time let pi, i > 4, be the result of substituting x — cos (9, y = sm6 in the terms 
of / of weighted degree i, where x has weight 2 and y has weight 1. Thus writing 
c = cos 9, s = sin 9 we have 

P4 = \k\c 2 + b 2 cs 2 + c 4 s 4 = q(c, s) (see Proposition EJ 
p 5 = bic 2 s + C3CS 3 + d 5 s 5 

and so on. Writing r = ro + r\t + r 2 t 2 + r^t 3 + . . . as before, where each 7^ is a 
function of 9 only, we find on substitution 

r Q = V-P4, 

4j>4ri + TqP5 = 0, from which we solve for r\, 

^P4TqT2 + &Pat\ + hp^r^ri +P6 r o = 0' from which we solve for r 2 , 

and so on. The assumption of an elliptic cusp of Gauss guarantees that p/4 is 

nonzero for all values of 9. 

The same ideas can be used in principle to parametrize the sections near 
a hyperbolic cusp of Gauss. It is not difficult to show that in this case (where 
b 2 > 2K1C4) the expression p<± above vanishes for exactly four values of 9 in the 
range < 9 < 2ir. In increasing order these are of the form < 6\ < 82 < 63 = 
2tt — 9 2 < 94 = 2n — 9\ < 2ir. In the ranges from 9% to 9 2 and from #3 to #4 the 
expression p^ is negative and in the other two ranges it is positive. Then we can 
use 

x = t 2 r 2 cos 9, y — trsin9, z = — t 4 in the first two ranges and 
x = t 2 r 2 cos 9, y — trsin9, z — t A in the other two ranges, 

determining r exactly as before, except that Tq = — when p 4 < 0. Needless 
to say this case is much more delicate than the elliptic case since the branches 
are not closed and the relationship between the value of t and the closeness of 
fit obtained by taking say four terms of the approximation to r will be difficult 
to determine. An example is shown in Figure 0| 
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4 Applications to pre-symmetry sets and symmetry sets 

Let 7 be a smooth parametrized plane curve. We seek first the 'pre-symmetry set' 
which is the set of parameter pairs (u±, u 2 ) such that there exists a circle tangent 
to 7 at 7(111), 7(^2)- The pre-symmetry set is symmetric about the 'diagonal' 
u\ = u 2 . A convenient way to find these pairs is to look for solutions of the 
equation given by a scalar product 

(71 " 72) ■ (Ti - T 2 ) = 0, (6) 

where ji = j(ui) and Tj is the unit oriented tangent to 7 at 7*. These solutions 
are precisely the pairs required, together with any pairs where Ti = T 2 : parallel 
oriented tangent pairs. For a closed convex curve the only such pairs arise from 
ui = U2 and these are easily identified as the diagonal. When a closed curve has 
inflexions other pairs with Ti = T 2 will be present. We need to beware of these 
when interpreting the pre-symmetry set. The pre-symmetry set contains much 
information about the symmetry set. Crossings of the diagonal correspond to 
endpoints of the symmetry set and horizontal or vertical tangents correspond to 
cusps. 

Secondly we want to plot the centres of the bitangent circles identified above. 
There are many ways to do this but here is a convenient one, couched in the 
language of complex numbers. Let c be the centre of the circle. Then, for some 
angle 6, 

(72-c) = e ie ( 7 i-c), T 3 = e M Ti. 

Then, as complex numbers, we have by dividing these equations c = (72 Ti — 
7iT 2 )/(Ti — T 2 ). Note that if Ti = T 2 , or, in practice, if these are sufficiently 
close together, and 71 ^ 72, then c is very far away and it will not appear on 
the diagram of the symmetry set. 

The above has been implemented in the Liverpool Surfaces Modelling Pack- 
age (LSMP, |Hj), also known as SingSurf. Note that it is necessary for 7 to be 
parametrized, which is the motivation for the theoretical results given in iJ21 We 
give below examples of umbilic points and elliptic cusps of Gauss (the latter 
obtained for the present by a different method). We shall give examples of hy- 
perbolic points elsewhere. See |4I5| for other methods and theoretical results 
related to the examples below. 

4.1 Umbilic points 

Elliptic points for which the principal curvatures as in are unequal make 
rather uninteresting examples since the intersection f(x, y) = k is a curve which 
is very nearly an ellipse, having exactly four vertices. The symmetry set has two 
smooth branches and the medial axis has one smooth branch, as for an ellipse. 
We concentrate here on the case of an umbilic point, where k\ — k 2 , and without 
loss of generality we can take these both equal to 1 by scaling the surface. Thus 
the equation of the surface has the form 



z = x 2 + y 2 + b x 3 + bix 2 y + b 2 xy 2 + b 3 y a + h.o.t. 
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Fig. 2. The umbilic case. Top left: several level sets f(x,y) — x 2 + y 2 + x 3 + 
2x 2 y + xy 2 — 3y 3 = k for k from 0.001 to 0.018. The curve becomes more circular 
as k decreases. The thick curves are the locus of vertices, obtained as in 0]. 
Notice that if k is too big then the level set ceases to be closed. Top right: The 
symmetry set and (boxed insert) the pre-symmetry set of the curve / = k for 
a small enough value of k (0.01 in fact) that the structure has stabilised. The 
heavily drawn part of the pre-symmetry set corresponds to the Y-shaped medial 
axis. The curve / = k was parametrized by the method of ^ using 10 terms. 
Bottom left: A closeup of the central part of the symmetry set; note the six cusps 
and two triple crossings. The medial axis is drawn heavily. Bottom right: The 
approximations to f(x, y) — 0.01 given by taking 1, 3, 4 and 10 terms as in 
the outermost curve is for 10 terms and is essentially indistinguishable from the 
exact solution. 
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By a technique explained in 0j we can find the locus of vertices of the family 
of curves f(x,y) = k for k > 0, and in there are results on triple intersections 
and cusps of the symmetry set. We shall not recall these in detail here but will 
point out how they are verified by the example. Clearly any rotation about the 
z-axis will leave the equation of the surface in the same form as above. If we 
rotate to make bo — b 2 then it can be shown that the three branches of the vertex 
locus make angles with the x-axis which are multiples of 60°, as in Figure top 
left. We pause here to establish the existence of such a rotation. 

Lemma 1. Let f be given by f(x, y) = x 2 + y 2 + box 3 + b\x 2 y + b 2 xy 2 + b^y 3 + 
h.o.t. Then we can rotate the coordinates in the x,y plane so that in the new 
coordinates u, v, the function f takes the form f = u 2 + v 2 + Bqu 3 + B\u 2 v + 
B uv 2 + B 3 v 3 + h.o.t. 

Proof. We rotate by an angle <j>. Expressing / in terms of the new coordinates 
(u,v), amounts to replacing x by u cos <f> + v sin cj> and y by — u sin (f> + v cos <f> in 
the expression for f(x,y). The new expression for / is of the form 
f = u 2 + v 2 + B u 3 + B\u 2 v + B 2 uv 2 + B 3 v 3 + h.o.t., where 
i?o = bo cos 3 (j) — b\ cos 2 (\> sin + b 2 cos 4> sin 2 ^ — 63 sin 3 4>, 
£>2 = 3bo cos 4> sin 2 <f> — b\ sin 3 </> + 2b\ cos 2 <f> sin <fi — 2b 2 cos <fi sin 2 <fi 
+b 2 cos 3 4> — 363 cos 2 4> sin (p. 
We need to show the existence of an angle </>o for which Bo = B 2 . Substitute 
U = tan(^) and write p = b 2 — bo, q — b 3 — b±. Then the equation B = B 2 
reads 

pU 6 + 6qu 5 - 15pU 4 - 20qU 3 + 15pu 2 + 6qU - p = 0. 

Note that the equation has pairs of roots of the form U,—l/U, corresponding 
to solutions (/>,(/> + n. Of course, if p = we can take = 0. Otherwise the 
discriminant of this degree 6 equation is a positive constant times (p 2 +q 2 ) 5 and 
hence > 0. For a degree six equation this means that there are either two or six 
real solutions for U. □ 

Remark We do not know a geometrical interpretation of the distinction between 
two and six solutions here. Note that these correspond to respectively one and 
three pairs of solutions </>, </> + 7r. There are a number of situations where through 
an umbilic pass one or three geometrically defined curves, for example ridges 
and sub-parabolic curves; see The present case appears to be different from 
these. 

In Figure H we show an example to illustrate the above methods. It confirms 
the results of 0] and for small enough k > 0, namely: (a) there are six vertices 
on the level set, resulting in three branches, of which one connects a maximum 
to maximum of curvature, one a minimum to minimum and one a maximum to 
minimum, (b) there are two triple crossings on the symmetry set, (c) there are 
six cusps on the symmetry set. Note that full detail is given on the symmetry 
set since the polar approximation is smooth and extremely close to the level set 
f = k. 
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Fig. 3. Elliptic cusps of Gauss, examples 1 (left) and 2 (right). The closed curves 
are the level sets f(x, y) = k, the other thin curves are the loci of vertices of level 
sets and the thick curves are the loci of inflexions. Notice that for small enough 
k there are four vertices (as predicted by @]) but as k increases the level set is 
first tangent to a branch of the locus of vertices not local to the origin and then 
intersects this locus in two points, giving six vertices. 

4.2 Cusps of Gauss 

We consider here elliptic cusps of Gauss, as in |J5J, and assume K\ > so that 
the level set f — k is, locally to the origin, a simple closed curve for small k > 0. 
According to the general results of ^ there are four vertices and two inflexions 
on this closed curve. (The situation for hyperbolic cusps of Gauss — and indeed 
for hyperbolic points — is much more complicated, with the pattern of vertices 
and inflexions giving rise to several cases, which are detailed in 4 .) Figure 01 
shows several level sets together with the loci of vertices and inflexions, for two 
elliptic cusps of Gauss. 

Example 1 is f(x, y) = x 2 + 2x 2 y + xy 2 + 2x 2 y 2 — xy 3 + y 4 , 
Example 2 is f(x, y) = x 2 + 2x 2 y + xy 2 + 2x 2 y 2 — xy 3 + &y A . 

Figure left, shows the curve of Example 1 together with approximations 
obtained as in In this figure the horizontal scale is exaggerated to improve 
clarity (otherwise the approximations and the actual curve are too hard to dis- 
tinguish!). We also show a hyperbolic cusp of Gauss example in Figure^] right, 
though we do not go on to examine the symmetry set here. 

Theoretical results on the symmetry sets of these curves appear to be much 
more difficult to obtain than those for elliptic, hyperbolic and ordinary parabolic 
points ( |4l5p . Thus evidence gathered from examples is all the more valuable. We 
show in Figure[3]the symmetry sets and pre-symmetry sets for the two examples 
of Figure |3J In both cases the stable situation as k — > is two smooth branches, 
one of which extends to infinity. Thus the structure of the symmetry set for 
an elliptic cusp of Gauss appears to be very much simpler than for an umbilic 



12 Diatta, Giblin, Guilfoyle and Klingenberg 




Fig. 4. Left: the thin curve is a level set f(x, y) = k for an elliptic cusp of Gauss 
as in Example 1 in the text, and the thick lines are obtained by taking r = tq 
and r = tq + nt in the approximations described in [|3J The approximation 
obtained by taking three terms r — r$ + nt + T2t 2 is visually indistinguishable 
from the true level set. Right: the thin curves are two level sets f(x,y) = k 
for values of k with opposite signs in the case of a hyperbolic cusp of Gauss 
f(x,y) = x 2 + ixy 2 + 3xy 3 + y A . The thick curve is obtained with the first 
term r = tq of the approximation of Clearly there are more problems here 
with capturing all the local features of the level sets. Note: in these figures the 
horizontal scale has been exaggerated compared to the vertical scale, as an aid 
to clarity. 

point. In particular there are no cusps and no triple crossings. Furthermore, the 
only distinction between the symmetry sets of sections of a surface close to the 
tangent plane at an ordinary elliptic point and at an elliptic cusp of Gauss is that 
in the latter case one of the branches of the symmetry set extends to infinity. On 
the other hand, the distinction between an ordinary elliptic point and an umbilic 
is very striking: in the latter case there are six cusps and two triple crossings as 
in Figure |21 

5 Conclusion 

We have given a simple method of parametrizing closed level set curves f(x, y) = 
k to any desired degree of accuracy and used it to compute symmetry sets and 
medial axes of sections of a surface z = f(x,y) in Monge form close to the 
tangent plane at an umbilic and an elliptic cusp of Gauss. The method applies 
also to level sets f(x,y) — k which are not closed; here the patterns of vertices 
and inflexions on the plane sections become much more complicated fl] and there 
are many different cases. Furthermore there are technical problems in ensuring 
that the approximations capture all the necessary local information about the 
level sets. We shall pursue these cases elsewhere. 

There are many related questions; for example, when a point on a surface 
moves from the hyperbolic region to the parabolic curve, how does the family 
of symmetry sets of the parallel plane sections behave? This is tantamount to 
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Fig. 5. Example 1 as in Figure [31 with annotated symmetry set and pre- 
symmetry set. The dashed section of the pre-symmetry set (boxed inserts) is 
due to the parallel tangents arising from the inflexions on f(x,y) — k; see 21 It 
is ignored when plotting the symmetry set (main diagram). The arrows indicate 
that a branch of the symmetry set goes to infinity when the bitangent circle 
becomes a bitangent line. The other endpoints are actual. Top: k is large enough 
for the level set to have six vertices. The additional branch of the pre-symmetry 
set is the loop towards the bottom left of the box and the corresponding piece 
of symmetry set is enlarged in the circle. As k decreases the number of vertices 
reaches its stable value of 4, as in the bottom figure. The medial axis for both 
values of k consists of the branch interior to the curve and the part of the other 
branch going to infinity to the left. 
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Fig. 6. Example 2 of Figure^ Here k is taken small enough to make the level set 
have its stable number of vertices, namely 4. The symmetry set has very much 
the same structure as Example 1, that is two simple branches one of which 
extends to infinity. 

considering a one-parameter family of surfaces, that is the curves will now belong 
to a fwo-parameter family. This again is the subject of further work. 
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